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§ ! Abstract 
(N 

>■ ■ An integrable structure behind Witten-Dijkgraaf-Verlinde-Verlinde (WDVV) equations is 
identified with reduction of a Riemann-Hilbert problem for a homogeneous GL{N,C) loop 
group. Reduction requires the dressing matrices to be fixed points of a loop group automor- 
phism of order two resulting in a sub-hierarchy of gl{N, C) hierarchy containing only odd 
symmetry flows. The model possesses Virasoro symmetry and imposing Virasoro constraints 
ensures homogeneity property of the Darboux-Egoroff structure. Dressing matrices of the 
reduced model provide solutions of the WDVV equations. 
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^ ■ 1. Introduction 

Massive topological field theories can be classified locally by the Darboux-Egoroff system of 
differential equations: 

O 

^ ' d ^ d 

^ • OUp; k=l k 

^. ■ A. = , (2) 

, , , where in addition one also assumes that 
<> 

^ d 

This was observed by B. Dubrovin in p[, 0. The jSij, 1 < i, j < N with i ^ j are called 
the "rotation coefficients". To be more precise, these Darboux-Egoroff equations are the 
compatibility conditions of the following linear system of differential equations depending on 
a spectral parameter A: 

d 

—Ajk{u,X) = f3ji{u)Aik{u,X) ; J = 1,...,N (4) 

N Q 

Y,-K-Ak{uA) = XAk{u,X) ] i = l,...,N (5) 
j=i '^'^j 

for each k = 1, . . ., A. Let V{u) = {(3ij{u){uj — Ui))i<ij<n and assuming that this matrix 
is diagonalizable, one can construct a local semisimple Frobenius manifold, which provides 
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an algebraic formulation of a massive topological field theory as follows. Find an invertible 
power series 

A{u, A) = Aoiu) + Ai{u)X + A2iu)X^ + • • • , (6) 
which is a solution of this system, then 

^ d 

Y^UkT^A^iu) = V{u)A^{u). (7) 

It is then straightforward to check that Ao{u)^^V{u)Aq{u) is a constant matrix, hence there 
exists an invertible complex matrix S such that 

N 

/i = J2^^^E^^ = S'' Aoiu)-'V (u) Aoiu) S (8) 
1=1 

is a diagonal matrix. Now let A{u) = Ao{u)S with A{u) = {aij{u))i<ij<N, then 

N Q 

Y^u,—A{u) = A{u)i^. (9) 

Dubrovin then showed that there exists on the domain Ui ^ Uj and 011021 ■ ■ -(Ini 7^ 
a local semisimple Frobenius manifold with scaling dimensions fia — yWi (see also 0). The 
diagonal metric is given by 

n 

d-s^ = Y,cili{u){duif. (10) 

j=i 

The flat coordinates are determined by 
where 

n 

Val3 = Y0iia{u)aip{u). (12) 
i=l 

Finally, the structure constants 

d^F{x{u)) ^ aia{u)aip{u)ai^{u) 
= dx^dx^dx^ = g ^) 

are the third order derivatives of the WDVV prepotential F{x). 

In this paper, we use a Riemann-Hilbert problem for the loop group of GL{N) to obtain 
a similar structure. The symmetry conditions (0) for the rotation coefficients are obtained by 
using the Cartan involution on the loop group. The model possesses Virasoro symmetry and 
imposing Virasoro constraints ensures the homogeneity property (^ of the Darboux-Egoroff 
structure. 
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2. Symmetry flows, Riemann-Hilbert problem and the r function 

We will introduce the G = GL{N) symmetry flows in terms of (extended) Riemann- 
Hilbert problem involving two subgroups of the Lie loop group G defined as : 



G^ = \ge G\g{\) = 1 + ' G+ = \ge G\g{X) = Y.g^'^ \ 



(14) 



here gi has grading i with respect to a homogeneous gradation defined by a derivation Xd/d\. 
It also holds that G+nG^ = I. Let the loop algebra corresponding to G be ^ = gl{N). This 
algebra splits into the direct sum G = Q+ ® G- with respect to the homogeneous gradation, 
where Q± are Lie algebras associated with the subgroups G±. 
Let the "un-dressed" wave (matrix) function be 

(Af oo \ N , , 

EE^]f-r =E^..e^-^""^ (15) 
j=l n=l J j=l 

We now define an extended Riemann-Hilbert factorization problem for the homogeneous 
gradation : 

^o(u, \)gi\) = e-\n, A) M(u, A) (16) 

where g : ^ G^G^ while G G^, M G G^. We use the multi-time notation with 
(u) = (ui,...,un) to denote multi-flows uj. Each argument uj, j = 1,...,A^ is an 
abbreviated notation for the multi- flows with n running from 1 to oo. 



For the symmetry flows one derives from (piq) expressions: 

d 



d 



^e(u,A) = -(04")0-^)_0(u,A) (17) 



du 



^M{n,X) = (0E];)0-i) M(u,A) (18) 



'j 



where (. . .)± denote the projections onto Q±. Note, that J2j=i dQ/duj = for any n > 0. 

In the following we will denote for brevity Uj = and dj = d/duf^ for j = 1, . . ., A^. 
It is natural to associate to the Riemann-Hilbert factorization a one-form : 

^ = - E E ResA ftr (Q-'^E^f)) duf^ = -Res, [tr (©"'^^o ^^^o) ) (19) 



j=l n=l 



which due to equation (|T^ is closed : dJ = with d = J2i=i J2'^=i duf^ Hence, locally, 

J can be written as of a zero-form being some function of the dynamical variables of the 
model. Conventionally, this function is denoted as the logarithm of the r-function: 



J = -d\o^- = -llt^-^ duf (20) 
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This equation can effectively be integrated to yield an explicit expression for the r-function. 
To do this it is required to extend the formalism from the loop groups to the centrally 
extended groups. As in case of the AKNS hierarchy with the r-function defined in we 
embed the framework in the centrally extended group with elements denoted as g and propose 

r(u) = (0|C^|0) (21) 

where |0) is the highest weight state such that X>o|0) = and (0|X<o = and furthermore 
(0|A;|0) = 1, where k denotes the direction associated to the central element within the affine 
Kac-Moody algebra. 

Let 9^~^^ be the grade —1 term of expansion: 

e = 1 + ^(-i)A-i + ^(-2U-2 + ... . (22) 

From equations (|l^) and (pOD (for n = 1) one finds : 

{9^-%, = -d,\ogT. (23) 

Introduce now the "rotation cefficients" Pij with i ^ j as the off-diagonal elements of 9^^^^ : 

P,^ = {9^-%^, tj^j = l,...,N (24) 

For 9^~^'' matrix we find from (|r^) : 

9/-^) = [E,,,^(-2)] +[9(-'Ke,,]9^-'^ (25) 

from which it follows that : 

0,(9^-%, = (9^~'^)..i0^-'^)jk, 3^1, 3 ^k. (26) 

Accordingly, the "rotation coefficients" [3ij given by (p^ ) satisfy (p. 

As found in the matrix elements of the dressing matrix = (%)i<i,j<Ar are given by: 

^.■(u,A)^4^x j (27) 



r u 



where the multi-flow arguments uj were shifted according to (uj — [A ^]) = {uj^^ — 1/A, u^'^^ 
1/2A2,...). 

The matrix coefficients of the wave matrix \1/: 

^(u,A) = 0(u,A)^o(u,A) (28) 

and the matrix coefficients of M(u, A) = \I^(u, A)5'(A) satisfy the linear system (§)-(|^), which 
are the compatibility conditions of equations (|l]). This is the Darboux-Egoroff system pro- 
vided we also impose the symmetry condition (Q) on the rotation coefficients. To do this we 
reduce the integrable hierarchy associated with the flow equation equations (|1^,(|18D using 
the automorphism a : 



a 



(X(A)) = ((X(-A))^)"\ XeG = GL{N) (29) 



which leaves the evolution equations (|T^-(0) invariant for odd flows only (labeled by n 
being an odd integer). Accordingly, we define the integrable sub-hierarchy by constraining 
the dressing matrices 0(u, A) and M(u, A) to the fixed points of the loop group automorphism 

e-^(u,A) = 0^(u,-A), M-i(u,A) = M^(u,-A) (30) 

with 9(u, A) and M(u, A) depending only on odd coordinates u: ( uf^~^^\ A; = 0, 1, . . .). The 
odd flows of the reduced sub-hierarchy preserve the conditions (^0|). The fixed points of the 
automorphism a form a subgroup of G = GL{N), called a twisted loop group of GL{N) 
(see e.g. 0). For the first term 9^~^^ of the expansion (|22|) the constraint (|30| ) implies that 
= 9^~^^'^. Accordingly, the "rotation coefficients" Pij associated to the subhierarchy 
via ( p4D are symmetric and give rise to the Darboux-Egoroff metrics. 

3. Condition of conformal invariance and homogeneity 

We propose here the Virasoro constraints in the setting of the Riemann-Hilbert problem 
approach to the integrable models. In this setting the Virasoro constraints are seen to arise 
as symmetries of the hierarchy. 

The action of the Virasoro symmetry flows on the dressing matrices is given by: 

5re(u,A) = ieL,e-')-e = iehe-')-e-j:T.nuf-j^ (31) 

j=l n=l OU) 

N 00 r) A/f 



6lM{n,X) = -(eL,e-i)+M = -(e4e-i)+Af-5:^n«f-— ^ (32) 

j=l n=l OU) 



where 



J N 00 

-^'■^':7r + EE^"^r^ k>-l (33) 



rl °° 

N 



= diag(4"V..,nSv"0=E«S"^^^^ 



i=l 



and Ik = —X d/dX. Both Ik and Lk for k > —1 generate a subalgebra of the Witt (centerless 
Virasoro) algebra: 

[lk, Ir] = {k- r)lk+r, [Lk, Lr] = {k - r)Lk+r (34) 

The signs in eqs. (|3l|) and (|32D have been chosen in such a way that the transformation 
also satisfies the Virasoro algebra [5^',^^] = {k — r)6^_^_J. when applied on the dressing 
matrices G and M. Furthermore, the transformations 5^ commute with the flows d/d- 



•u 



(n). 



[6l-^]Q{u,X) = (35) 



d 



du 



(n) 
3 



Accordingly, the Virasoro transformations from eqs. (|3ll) and ( ^2t ) constitute additional 
symmetries of the integrable model defined by the Riemann-Hilbert factorization problem. 
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One can show that only the even Virasoro flows will preserve the conditions (^) and define 
additional symmetries of the reduced sub-hierarchy associated with the twisted loop group 
of GL{N). We can therefore consistently impose the constraints (5^0(u, A) = with even 
k. For k = the equation (|3l| ) specializes to: 



S^e{n,X) = X^e-j:T.nuf^ (36) 

as follows from identity — (GA(iG^^/(iA)_0 = XdQ/dX. Hence for k = the constraint 
(5^ 0(u, A) = takes a form 

(^A^-Eo)e(u,A) = (37) 
where we have introduced the Euler vector field Eq: 

N oo o 

^o = EE-^rAT (38) 



j=l n=l du 
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The homogeneity condition (|37|) is also a condition for the conformal Frobenius manifold. 
By projecting in equation (0) on the terms with —1 grade we obtain : 

5o^^(-) = -^(-)-EE-r^ (39) 

j=ln=l OUj 

where 9^~^^ is the grade —1 term of expansion: (P^. 
Hence, the condition (^) implies that 

Eq ^(-1) = -^(-1) (40) 



Plugging relations (p3|) and (El) into equation (p9D we obtain 



^ °° , glog r 
j=i „=i 



^o^a.logr = 9.5o'^logr = -9.EE^"^^^' z = l,...,iV (41) 



N oo fiO 

Oh = -A/-EE^"^^ (42) 



j=i n=i 



4. Virasoro flows on the r function 

This subsection studies action of the general Virasoro transformation on the r-function of 
the underlying integrable system. 

Recall an expression for derivatives of logr given in ([T9|) and (pOf). Applying di on both 
sides of eq. ( pOD we obtain an identity: 

-^^ = Resx{tT{eE^fe-'Eu)) i,j = l,...,N, n>Q (43) 



duj 'du. 
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Combining eqs. (|3TD and (p!^)-(]20|) we find 

sf-^ = Res. (tr (e4»'e-A(ei,e-i)_JJ (44) 



Since \d{\ ^)/dX = —A ^ the last equation becomes 

N 00 

Sld.logT = -ResA (tr (^,,6^^9-1)) "EE nu^Res^ (tr (eEg+^^e-^E,,)) (45) 

j=l n=l 

Using the identity (^3D we rewrite the last equation as 

TV 00 o 1 

5ld,\ogT = -Res. (tr (E^.e/.e"^)) -EE nufd,^^ (46) 

j=l n=l OUj 



The cases = 0, 1 were calculated and eq. (|45|) gives : 

,{«) 



iV 00 o 1 

a.5riogr = -9.EE^"^|^^ ^ = 0,1 (47) 



a.5r,iogr = -5.EE^"^|^- (4^ 



j=l n=l 

and 

Af 00 
j=l n=2 <9m^ 

Only k = case applies to the reduced case. 
From eqs. (|3T| ) and (^2]) we find : 

sl(e-^M) = -Q~\eLke-^)_M -Q-\QLke-^)_M = -Lk(e-^M 

J N 00 

= A'=+^;^ (q-'m) -EE ^"^^ir'^ (Q-'m) (49) 

•^-^ j=l n=l 

Recalling the factorization problem (0) we find 

<5r(e-^M)=v]/o(u,A)A^+^^ (50) 

since 

J N 00 

X'^'^^oin, A) = E E ^"^4"^'^^o(u, A) (51) 

"'^ i=l n=l 



Recall, the definition (pTD of the r-function Tg (u) = (0|\E'o5' |0), rewritten here in the form 
which explicitly shows dependence of r on the loop group element g{X)- 
It now appears that 



S^rg = m'^'^o^ |0) = r,.+i,,/,A (52) 
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So under infinitesimal diffeomorpism g g + X^^^dg / dX the r-function transforms as 

Tg Tg+xk+Ug/dX =Tg + Tg (53) 

5. Virasoro symmetry and dressing 

Let 

Sf = -^-E^\ 6, = 6f = d, - XE,,, j = l,...,iV (54) 

OUj ' 

be a set of commuting operators : 

[6f , 6l'^] = 0, (55) 
The above operators and the Virasoro operators Lk from (|3^) commute with each other: 

[5f\Lfc] = 0, z,j = l,...,N (56) 
Both operators ( [5^ ) and ( |55D annihilate the bare wave function from (|TB|) : 

Lfc$o = 0, (5f^^o = 0, j = l,...,iV. (57) 

Similarly the operators: 



satisfy the commutation relations (^) and (Q) and annihilate the wave (matrix) function 
^(u,A) defined in (^ : 

Ck^ = 0, VP^ = 0, j = l,...,N. (59) 



Since 



eSj-e-i = dj + eidjQ-^) = dj + (eA%e-i)_ (eo) 

we obtain 

V, = vf = e {d, - XE,,) = d, - {QXE,,e-') ^ = a, - a% - v, (ei) 

where 

V, ^ \d^-'\ E,j], iV,)ki = i6,j-6k,)/3ki (62) 

for ^(-1) defined in (H). 
Consider 

Recall, now from (|36D that the condition (5]Lq9(u, A) = (|37D implies 

i=ln=l 
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and accordingly 

A^ = (0EE^"^^]f0'V^ (65) 

"-^ j=l n=l 

Let us now study the operator Ck=~i = 6L_i0^^ from (^) for satisfying the homo- 
geneity condition 



dX^~dX 



= -i^ + ^i^e-' + EEnur^E^r^-^ (66) 

j=l n=l 



^ N oc N c<y 

dX 



jr=l n=l j=l n=l 



J N oo 

.J- + A-EE-r(ei^]fe-- 

"'^ j=l n=l 



and for = 0, n > 1 one obtains 



£_i = + f/ + \-'[0^-'^ ^u] = -^ + U + A-V (67) 



The components of the matrix V = [9^ U] are 



Vij = {uj - Ui)elj = {uj - Ui)f3ij, i,j = 1,...,N (68) 

It follows from the construction that the wave (matrix) function \E'(u, A) defined in (^) is 
annihilated by £_i and Vj 

C-i^iu, A) = ^ = {U + A-V)* (69) 
uA 

and 

P,*(u, A) = ^ 1^ = (A% + Vj)^iJ . (70) 

Compatibility of the above equations amounts to [Dj , C-i] = 0, which follows by dressing 
of [6j , = 0. From commutation relations [Dj , C-i] = and [Di , Pj] = one obtains : 

d,V = [V,,V] (71) 
djV, = d,V, + [V,,Vi\ (72) 
[V,E,,] = [V,,U] (73) 

6. Isomonodromic Tau function ri 

We now consider the coefficients —dj log r of the one-form J' from (|19]) with the conformal 
constraint (|64D imposed. Recall, that this constraint follows from the homogeneity condition 
(p7|) and in view of relations (^Tj) and ( ^21) it holds in the reduced case that Eqt = fir and 
EoPki = —Pki where Eq = YljLi Ujd/duj and /i is independent of variables Uj,j = 1, . . ., N. 
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Plugging relation (|6^ into the defining equations (|T9D and ( pOD we obtain for n = 1: 

9,logr = -ResA |^tr (^Ej.iej^flnu^I'^E^fQ'')-^! , J = l,...,N (74) 

(k) 

In the rest of the paper we set all u] = for A; > 1 leaving the model defined only in terms 
of the canonical coordinates Uj,j = 1, . . .,N and assume from now on that a{g{X)) = g{X) 
and that {^) holds. 

The formula ( |7^ simplifies significantly and takes the form of 

9,logr = -Res, (tr (Ej.iQXUQ-')^)) = -tr (e,, ([e^-'\U] + ^[e^-'\ [e^-'\U]]^^ 

(75) 

Using invariance of the trace tT{A[B, C]) = tr([74, B]C) we obtain expression: 

dj\ogT = ^tT{V,V) (76) 

with Vj, V matrices defined in (^) and ([68| ) and related through the commutation relation 

The equation ( ffTD can be cast into the Hamiltonian form (see also [Q) 

d,V = {V,H,} (77) 
with respect to the standard Poisson bracket {■ , ■} on so{N) and with the Hamiltonians 

H^ = tJ2^^_ (78) 

The commutation relation ( [73D in components gives a relation 

{uk-Uj){Vi)kj = {5ik- 5ij)Vkj, k^j = l,...,N (79) 
which can be used to prove that 

l tr im = E = - E = _2H, (80) 

From equations (|7lD and ( [72l) it follows that dtr X^ili (^^) c^^^i = and the one-form 
J2iLi Hidui is closed. Locally, this one-form can be represented by the derivative of the 
logarithm of the so-called isomonodromy tau-function r/ 

N ^ 

dlogTi = ^Hidui = --dlogr. 
1=1 ^ 

which implies that upto a multiplicative constant : 

ri = 1/VT. (81) 
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This formula was also obtained in a different way in 0. 
7. Frobenius structure 

Expand, the dressing matrix M{u, A) according to the grading : 

M{u, A) = Mo(m) + Mi(m)A + M2(m)A^ + . . . (82) 

Then, according to the condition (pO|) the zero grade term Mo(m) = M{u, 0) = {mij{u))i<ij<N 
is an orthogonal matrix Mq{u) = Mq'^(u) which as seen from eqs. (|TB|) and ( P^D also satisfies 
the flow equations : 

djMoiu) = [e^-'\u) , Ej,]Mo{u) = Vj{u)Mo{u), S^Mo{u) = ~V{u)Mo{u) (83) 
and hence also 

N 

Y,u,d,Mo{u) = V{u)Mo{u). (84) 

Define 

V = Mq-VMo, (85) 

it is easy to see that V satisfies: 

a,V = 0, V^ = -V, CV = 0, j = l,...,N (86) 

and accordingly V is a constant. We can bring V into Jordan normal form fi = fi + Rq, with 
jl its semisimple part and Ro its nilpotent part, i.e. there exists a complex invertible matrix 
S, such that 

V = SfiS~\ (87) 

Now define 

M(m) = Mo{u)S = {m,,{u))i<i,j<N, (88) 

then using ( ^31) we obtain: 

djM{u) = [e^-^\u) , ^jj]M(m) = Vj{u)M{u), 5^M{u) = -M{u)n (89) 
and hence also 

TV 

UjdjM{u) = M{u)fx. (90) 
Now assume that V is diagonalizable, i.e. the matrix V{u) has no nilpotent part, then 

N 

fJ' = 'YfiiEii, with fii = -fiN+i-i (91) 
1=1 

and the columns of M{u) are eigenvectors for the Euler vectorfield Eo. Hence, on the domain 
Ui 7^ Uj and miim2i ■ ■ ■ mTvi 7^ we have a local semisimple Frobenius manifold (see e.g. 0]) 
with Lame coefficients 

h, = m,i. (92) 
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Define, 

V = Mi<a,f3<N = M^M = S'^S, and denote = {t%<a,p<N, (93) 

then fir] + rjfi = and the derivatives with respect to the fiat coordinates x", a = 1, . . .,N 
are given by 

and 

E ^"Vim.,--. (95) 



' o,p=l 



The metric equals 



N N N 

ds^ = '^h1{duif = E Vapdx'^dx^, with r^gg = mjamig. (96) 

i=l a,P=l 1=1 

and the structure constants are given by 

AT 

^ I 'He 



~ ■ (97) 



Since we have constructed everything in such a way that the columns of M{u) are eigenvec- 
tors, our Frobenius manifold has scaling dimensions — f^i- 

Now, multiply the wave (matrix) function \E'(m, A), defined in (p8|), from the left by 

M-\u) = S-^M^\u) = r]-^S^M^{u) = r]-^M^{u) (98) 

and from the right by S to obtain the wave function y{u, A): 

y{u, A) = M-\u)^{u, X)S = S-H'Io\u)M{u, X)g'\X)S . (99) 

y{u, A) satisfies, as in ^ : 

y{u,-xfr]y{u,x) = T] (100) 
dy{u,x) 



dui 

dy{u,x) 

dx"' 

dy{uA) 

dX 



xs,{u)y{u,x) (101) 

XC^{u)y{u,X) (102) 
{U{u) + A- V) y{u, X) (103) 



where 



S,{u) = M-\u)EuM{u) (104) 
U{u) = M-\u)UM{u) (105) 
/i = M-\u)V{u)M{u) = S-^VS (106) 
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and 

N N 

Caiu)= Y.^lpEp,= E (107) 

/3,7=1 /3,7,5=1 

For general V{u), which may contain a nilpotent part, we have ^^i] + r]fi = 0. Define 
again y{u, A) by (pOj) then also 



Miu, A) = M-\u)M{u, X)S = I + 0(A) (108) 
such that y{u, A) = M{u, X)S-^g-\X)S and 

Miu,-\friM{u,\) = r]. (109) 
From equation (|103|) one obtains : 

g-^^^^SM-H^-UM-^^M)S-^ (110) 

Assume, that g = g-g+ then 

Define G = g+SM'^S-^ then 

(71^^ = -jGoSfiS-'G,' = -^GoVGo^ (112) 

with matrix Gq = {g+S^A^^S^^)o = {g+)o being orthogonal for a fixed point g of the 
automorphism a. 
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